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SECTION  I 

PROGRAM  DESCRIPTION 

1 . 1  INTRODUCTION 

The  Guidance  and  Control  Systems  Division  of  Litton  Systems,  Inc. 
is  involved  in  a  multi-year  program  in  the  research  and  develop¬ 
ment  of  a  nuclear  magnetic  resonance  (NMR)  gyro  for  use  in  iner¬ 
tial  navigation  systems. 

In  the  magnetic  resonance  gyro,  angular  craft  rate  is  sensed 
through  a  shift  in  the  precession  frequency  of  an  ensemble  of 
atomic  or  nuclear  magnetic  moments.  The  Litton  NMR  gyro 
utilizes  ensembles  of  two  noble  gas  isotopes  in  order  to  elimi¬ 
nate  explicit  dependence  on  the  magnetic  field.  In  order  to 
obtain  the  gyroscopic  information  in  this  scheme,  the  noble  gas 
ensembles  must  first  be  polarized  so  that  coherent  precession 
about  an  applied  magnetic  field  can  take  place  and,  secondly,  a 
means  for  detecting  the  nuclear  precession  must  be  provided.  in 
the  Litton  approach,  the  noble  gas  ensembles  attain  a  net  mag¬ 
netic  moment  through  spin  exchange  with  an  optically  oriented 
alkali  vapor.  In  addition,  the  alkali  vapor,  through  a  magne¬ 
tometer  effect,  detects  the  weak  magnetic  fields,  which  are 
generated  by  the  precessing  nuclear  moments  of  the  two  noble 
gases,  and  thus  provides  the  gyro  angular  rate  information. 

The  purpose  of  the  physics  research  program  is  to  develop  a 
stronger  theoretical  and  empirical  understanding  of  the  alkali- 
atom  noble-gas-atom  system.  Of  particular  interest  are  the  spin 
exchange  interactions  between  the  alkali  valence  electrons  and 
the  noble  gas  nuclei,  the  noble  gas  nuclear  spin  relaxation 
mechanisms  as  they  pertain  to  the  NMR  gyro  and  other  related 
phenomena  which  impact  the  NMR  gyro  technology. 
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1.2  OBJECTIVES  AND  STATUS  OF  THE  CURRENT  PROGRAM 

The  research  tasks  described  in  Sections  II  and  III  of  this 
report  are  part  of  our  continuing  effort  to  understand  and  param¬ 
eterize  the  alkali-atom  noble-gas-nucleus  spin  exchange  inter¬ 
action.  In  addition,  investigations  of  the  spin  exchange  process 
also  provide  opportunities  to  study  the  pertinent  relaxation 
mechanisms  of  the  noble  gas  nuclei. 

Reported  in  Section  II  are  the  results  of  the  theoretical  and 

experimental  investigation  on  the  relaxation  of  polarized  noble 

gas  nuclei  through  the  interaction  of  the  nuclear  quadrupole 

moment  with  the  electric  field  gradient  arising  from  collisions 

with  the  walls.  In  the  last  AFOSR  annual  Technical  Report,  a 

simple  theoretical  model  of  this  effect  was  developed  and  was  in 

agreement  with  the  experimental  data  only  qualitatively.  No 

quantitative  analysis  of  the  data  was  possible  with  the  model. 

In  the  present  work,  the  model  has  been  extended  to  incorporate 

the  second  order  expansion  of  the  nuclear  density  matrix.  Under 

131 

this  new  model,  measurements  of  the  relaxation  of  Xe  nuclei 
were  analyzed  quantitatively,  and,  as  a  result,  the  sticking  time 
of  the  Xe  atom  onto  the  container  wall  and  the  strength  of 
the  nuclear  quadrupole  interaction  have  been  estimated  for  the 
first  time. 

In  Section  III,  we  report  on  our  measurement  of  the  shift  in  the 

Larmor  precession  frequency  of  the  noble  gas  nuclei  during  spin 

exchange  with  an  optically  oriented  alkali  vapor.  The  results 

confirm  that,  for  a  given  alkali  polarization  environment,  the 

difference  in  frequency  shifts  for  different  isotopes  of  the 

129  131 

same  noble  gas  atom,  such  as  Xe  and  Xe  ,  is  smaller  by,  at 
least,  an  order  of  magnitude  than  the  difference  in  frequency 
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129  83 

shifts  for  different  noble  gas  species,  such  as  Xe  and  Kr 
Experimental  results  are  presented  in  terms  of  the  parameters 
relevant  to  the  NMR  gyro. 

Investigations  into  this  area  of  physics  have  relevance  not 
only  to  NMR  gyro  research  but  also  to  other  possible  areas  of 
research  such  as  magnetometers,  atomic  clocks,  masers,  laser 
stabilization,  and  polarized  nuclear  targets.  It  was  in  this 
light  that  the  result  of  the  work  reported  in  Section  II  has 
been  prepared  as  a  scientific  paper  to  be  submitted  for  publica¬ 
tion  to  The  Physical  Review.  The  preprint  of  this  paper  is 
included  as  paragraph  2.2. 

In  the  course  of  this  work,  an  idea  was  conceived  to  coat  the 
container  wall  with  RbH  to  reduce  the  nuclear  quadrupole  inter¬ 
action.  Successful  implementation  of  the  idea  has  been  demon¬ 
strated  in  the  data  reported  in  Section  II.  We  believe  that  the 
invention  of  the  coating  material  and  coating  technique  is  a 
significant  step  forward  in  the  development  of  NMR  devices.  It 
is  for  this  reason  that  we  have  submitted  a  written  disclosure 
of  invention  to  the  company  for  consideration  in  applying  for  a 
patent.  The  patent  application  is  under  consideration  for 
further  action  by  the  corporate  headquarters. 

1.3  PROFESSIONAL  PERSONNEL  ASSOCIATED  WITH  THE  PROGRAM 

The  principal  investigator  of  these  research  efforts  has  been 
Dr.  Tae  M.  Kwon.  This  work  has  also  been  strongly  supported  by 
Dr.  Charles  H.  Volk,  and  Dr.  John  G.  Mark.  We  also  wish  to 
acknowledge  the  contribution  of  Mr.  Howard  E.  Williams  for  the 
design  and  Mr.  William  P.  Debley  for  the  construction  of  various 
components  of  the  apparatus,  including  the  resonance  cells  and 
lamps  used  in  this  research  effort.  Finally  we  acknowledge 
Dr.  Edward  Kanegsberg  for  many  helpful  discussions. 
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SECTION  II 

QUADRUPOLAR  NUCLEAR  SPIN  RELAXATION 

2 . 1  INTRODUCTION 

1  pq  ni  1  oq 

Two  nuclear  species,  Xe  (1=1/2)  and  Xe  (1=3/2),  or  Xe 
8  3 

and  Kr  (1=9/2)  are  required  under  Litton 's  approach  to  the 

development  of  a  nuclear  magnetic  resonance  (NMR)  gyro.  Asso- 

83  131 

ciated  with  the  nuclei  of  Kr  and  Xe  ,  and  characteristic  of 
nuclear  spins  greater  than  1/2,  is  the  nuclear  quadrupole 
moment  which  interacts  with  electric  field  gradients.  Field 
gradients  arise  at  the  nuclear  site  during  collisions  of  noble 

gas  atoms  with  the  walls  of  their  containers.  This  interaction 

131  83 

is  more  pronounced  in  the  Xe  nuclear  ensemble  than  in  Kr 

Attention  is  directed  to  the  fact  that,  limited  by  this  inter- 

131 

action,  the  transverse  relaxation  time  (T2)  of  the  Xe  nuclei 
published  in  available  literature  prior  to  this  contract  had 
been  on  the  order  of  a  few  hundred  milliseconds.*  At  the 
initial  stage  of  this  contract,  it  was  improved  to  5  to  10 
seconds  in  15  ml  spherical  samples.  The  extremely  short  T2  made 
it  difficult  to  improve  our  understanding  on  the  subject  of  the 
alkali-atom-noble-gas-nucleus  spin  exchange  interaction,  which 
is  essential  for  the  development  of  an  NMR  gyro. 

Both  theoretical  and  experimental  understanding  of  the  nuclear 
quadrupole  interaction  continues  to  improve.  In  this  section, 
we  report  T2  of  a  Xe  nuclear  ensemble  as  long  as  220  seconds. 
This  improvement  was  accomplished  by  coating  the  cell  wall  with 
RbH .  With  the  experimental  advantage  of  a  larger  T2 ,  the 
parameters  relevant  to  the  nuclear  quadrupole  interaction  were 
estimated  for  the  first  time.  The  results  of  investigations 

*D.  Brinkham,  E.  Brun,  and  H.  H.  Staub,  Helv.  Phys.  Acta  3J3  431 
(1962) 
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have  been  prepared  as  a  scientific  paper  for  consideration  of 
publication  in  the  Physical  Review.  The  preprint  of  this  paper 
is  included  in  the  Appendix  to  this  section  in  lieu  of  a  para¬ 
phrased  report  on  the  research  performed  under  this  contract. 

The  developed  process  of  data  analysis  can  now  be  utilized 
routinely  in  the  gyro  development  to  determine  the  characteris- 
tics  of  cells  containing  Xe  (or  Kr  )  nuclei. 

Investigations  of  noble  gas  nuclear  quadrupole  moments  inter¬ 
acting  with  the  walls  of  their  container  will  continue,  particu¬ 
larly  for  RbH  coated  glass  walls.  In  order  to  improve  the 
uniformity  of  the  coating,  the  reaction  will  be  studied  in  glass 
containers  made  of  other  than  pyrex  glass.  Non-pyrex  glass  will 
be  studied  because  the  rubidium  reaction  with  hydrogen  molecules 
has  been  shown  to  be  sensitive  to  the  presence  of  impurities 
(such  as  water,  oxygen,  etc.)  normally  embedded  on  pyrex  glass.* 

2.2  PUBLICATION  PREPRINT 

A  copy  of  a  preprint  to  be  submitted  to  the  Physical  Review  for 
consideration  of  publication  follows  herein. 


*A.  Tam,  G.  Moe ,  and  W.  Happer,  Phys.  Rev.  Lett.  3j>  1630  (1975) 


2-2 


404577 


131 

Quadrupole  nuclear  s;£Ln  relaxation  of  Xe  in  the  presence 
of  rubidium  vapor 

T.  M.  Kwon,  J.  G.  Mark  and  C.  H.  Volk^ 

Litton  Guidance  &  Control  Systems 

5500  Canoga  Avenue,  Woodland  Hills,  California,  91365 
(Received 

Nuclear  spin  relaxation  of  '*'~5'1'Xe  produced  by  collisions  of 

the  xenon  atoms  with  the  walls  of  the  experimental  cell  has  been 

studied.  The  xenon  nuclear  ensemble  is  polarized  through  spin 

exchange  collisions  with  an  optically  oriented  rubidium  vapor. 

131 

Depolarization  of  the  Xe  ensemble  in  the  wall  collisions  is 

131 

assumed  to  be  due  to  the  quadrupole  coupling  between  the  Xe 
nucleus  and  the  electric  field  gradients  which  appear  at  the 
nucleus  during  collisions  of  the  atoms  with  the  cell  walls.  Model 
equations  describing  the  nuclear  relaxation  are  derived  and 
estimates  of  the  electric  field  gradients  are  given. 


PACS  numbers:  76.60.Gv,  76.60.Es,  07.58.+g 
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I .  INTRODUCTION 


We  recently  reported  on  measurements  of  the  transverse 

8  3 

nuclear  spin  relaxation  of  Kr  in  the  presence  of  rubidium 

vapor. *  The  decay  of  the  krypton  nuclear  spin  was  found  to  be 

dependent  on  the  orientation  of  the  experimental  cell  with 

respect  to  an  axis  defined  by  the  static  magnetic  field  present 

in  the  apparatus.  The  angle  between  the  cell  and  the  magnetic 

field  which  yielded  the  longest  decay  time  was  termed  the 

201  2 

'magic'  angle  after  similar  work  done  on  Hg  .  Furthermore, 
we  observed  that  the  decay  signal,  off  the  magic  angle,  possessed 
a  non-exponential  character.  These  effects  were  modeled  as  an 
interaction  between  the  nuclear  quadrupole  moment  and  the 
electric  field  gradients  present  at  the  krypton  nucleus  during 
collisions  of  the  krypton  atoms  with  the  walls  of  the  experi¬ 
mental  cell.  The  model  describes  the  transverse  nuclear  relaxa¬ 
tion  by  21  precessional  frequencies,  where  I  is  the  nuclear 

8  3 

spin;  I  =  9/2,  for  Kr.  We  are  able  to  show  a  qualitative 
agreement  between  this  model  and  the  experimental  measurements. 

A  quantitative  analysis  of  the  krypton  data  was  not  possible  due 
to  the  complexity  of  the  signal. 

We  report  here  experimental  data  on  the  transverse  nuclear 
131 

spin  decay  of  Xe  (nuclear  spin,  I  =  3/2) ,  which  exhibits 
relaxation  through  an  electric  quadrupolar  interaction  in 
collisions  with  the  cell  walls.  In  order  to  simplify  the 
experimental  data  analysis  we  have  re-formulated  our  model  of 
the  nuclear  quadrupolar  interaction,  viewing  the  precession  of 
the  nuclear  ensemble  as  21  oscillators  each  with  a  single 
frequency.  In  this  way  we  treat  21  two-state  Bloch  type 
systems^,  instead  of  treating  one  complex  system  of  21  frequen¬ 
cies.  This  procedure  is  similar  to  the  'ficticious  spin-1 /2 
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operators'  presented  by  Vega  and  Pines.  Furthermore,  our 
previous  model  is  extended  in  this  work  by  incorporating  results 
of  a  second  order  expansion  of  the  nuclear  density  matrix  as 
given  by  Cohen-Tannoud j i  for  the  case  of  relaxation  of  nuclear 
spin  3/2  under  a  quadrupolar  interaction.5  In  this  new  model, 
we  show  that  the  second  order  expansion  leads  to  a  coupling 
among  the  21  oscillators. 

II.  THEORETICAL  MODEL 


We  consider  in  this  section  the  transverse  relaxation  of  a 
nuclear  spin  ensemble  of  nuclear  spin  I  =  3/2  under  a  quadrupolar 
interaction.  We  will  rely  on  results  presented  in  Refs.  (1  and  5) 
and  repeat  only  selected  details  here,  as  clarity  demands. 


A.  The  energy  eigenvalues.  The  full  spin  Hamiltonian  of  the 
system  is  taken  to  be: 


H  =  Hm  +  hq 

where  HM  is  the  static  magnetic  interaction 


Hu  =  -  fiH  I_ 
M  o  Z 


(1) 

(2) 


with  the  nuclear  gyromagnetic  ratio,  Hq  the  external  magnetic 

field,  taken  to  be  along  the  z-axis,  and  I  ,  the  z-component  of 

z 

nuclear  spin.  Hq,  the  quadrupolar  Hamiltonian,  represents  the 
interaction  of  the  nuclear  electric  quadrupole  moment  with  an 
electric  field  gradient.  As  was  argued  in  Ref.  (1),  the 
electric  field  gradient  arises  during  collisions  of  the  noble 
gas  atoms  with  the  walls  of  the  experimental  cell.  If  we 
assume  that  for  a  given  collision,  the  resultant  electric  field 
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gradient  possesses  cylindrical  symmetry,  then  the  quadrupole 
Hamiltonian  can  be  written 

hq  -  4lm?rr  jl  13  °°s2e-»  K2  -  ki-d] 

Jz]  <3> 


3  2  +  2  _2 
+  |  sin  G  (I  +1  ) 

2  2 

where  e  is  the  electronic  charge,  eq  =  9  V/9z  ,  is  the  electric 
field  gradient  under  the  assumption  of  cylindrical  symmetry, 

C  is  the  quadrupole  moment  of  the  noble  gas  nucleus,  ©is  the 
angle  between  the  quantization  axis  and  the  cylindrical  axis  of 
the  electric  field  gradient,  and  1~  are  the  raising  and  lowering 
operators  associated  with  the  spin  operator,  I,  where  in  the 
cartesian  reference:^ 


+  2  sin©  cos© 


VI++I"> 


(I++l“) 


I1  =  ^  ±  ily  • 

As  in  Ref.  (1),  we  make  the  following  definitions: 


a  = 

(3cos2©  -1) , 

(5a) 

b  = 

sin©  cos  © , 

(5b) 

c  = 

sm  ©, 

(5c) 

a  = 

e2Qq/4I (21-1) , 

(5d) 

u>  = 

-7hH  . 
o 

(5e) 
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The  quadrupole  Hamiltonian,  H^,  is  small  compared  to  the 
magnetic  Hamiltonian  HM,  and  thus  we  compute  the  energy  eigen¬ 
values  for  the  full  Hamiltonian  employing  perturbation  theory®. 
The  energy  levels  for  a  nucleus  with  spin  3/2  are  computed 
through  second  order: 


E 


E 


E 


E 


1 


2 


3 


4 


3  3  27a2b2  ,  27a2c2 

2W  +  2aa - 5—  +  — 8^“ 


2"  “  2aa  ‘ 


2  2  2  2 
21a  27a  c 


oj 


8  co 


1..  3  .  27ot2b2  27cy2c2 
■2"-2aa  +  ~ « - OZT 


3 , ,  .  3  27a2b2  27a2c2 

2  W  +  2aa - 5 - SZT 


(6a) 


(6b) 


(6c) 


(6d) 


B.  The  system  equations.  The  evolution  of  the  observables  of 

the  system  are  found  by  solving  the  density  matrix  equation  of 

q 

the  ensemble,  which  is  given  by: 

If  -  K  Ip-  "I  I’l 

and  then  computing  the  expectation  value  for  a  system  observable, 
as : 


<0>  =  Tr  [p,  ol 


(8) 
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where  H  is  the  Hamiltonian  of  the  system  given  in  Eq.  (1) ,  p  is 
the  density  matrix  of  the  nuclear  ensemble  and  0  is  an  operator 
corresponding  to  an  observable  of  the  system.  For  the  situation 
represented  by  Eq.  (1),  in  which  there  is  a  large  time- 
independent  interaction,  H„,  and  a  much  smaller  but  time- 
dependent  term,  H^ft),  it  is  convenient  to  transform  the  density 
matrix  into  the  interaction  representation,  which  is  defined 

K...  10,11 


-  (i/h) HMt  (i/h)HMt 

Pit)  =  e  H  p*(t)  e  " 


where  p*  (t)  represents  the  density  matrix  in  the  interaction 
representation.  The  differential  equation  for  p*  is  easily 
found  to  be: 


ft 

1 k  *  K  HQltn 


where  H*  is  given  by: 


(i/h) H„t  - (i/h) HMt 

Hq (t )  =  e  "  HQ(t)  e 


The  time  dependence  of  Hq  is  taken  to  be  that  due  to  stochastic 
collisions  of  the  xenon  atoms  with  the  walls  of  the  experimental 
cell,  in  the  same  manner  as  that  treated  in  Ref.  (5). 
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In  the  limit  that  the  quadrupole  interaction  can  be  taken 

1 2 

as  a  perturbation,  we  solve  Eq.  (10)  in  an  iterative  manner  , 
which  through  second  order  is  given  by: 


-  K  Vu] 


\  2/"t  j 

r  *  * 

i  *  ) 

a ! 

|  p  (o)  ,  Hq  ( t ' ) 

].  Hq  ( t  )  j  df 

• 

(12) 

The  first  term  in  Eq.  (12)  was  shown  previously1  to  lead  to  a 
transverse  decay  of  the  nuclear  ensemble  through  a  dephasing  of 
the  precessing  Zeeman  levels.  The  second  term  in  Eq.  (12)  has 

5 

been  treated  by  Cohen-Tannoud ji  for  the  case  of  relaxation  of 
201 

Hg  under  a  quadrupole  interaction.  Since  the  nuclear  spin 
201 

of  Hg  is  3/2,  we  can  directly  apply  the  results  of  Cohen- 
Tannoud  ji's  work  to  our  analysis  of  1;11Xe. 

The  commutators  in  Eq .  (12)  are  evaluated  under  the 

following  assumptions: 

1.  The  first  order  commutator  can  be  evaluated  with  a 
diagonal  Hamiltonian  whose  eigenvalues  are  computed  through 
second  order;  and, 

2.  The  square  of  the  product  of  the  correlation  time  of 
the  interaction,  and  the  energy  level  separations  of  the 
nuclear  Zeeman  states,  is  much  less  than  unity. 

The  first  term  in  Eq.  (12)  is  then  evaluated  by  employing 
some  straightforward  matrix  algebra,  while  the  evaluation  of  the 
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second  terms  is  found  in  Ref.  (5).  Eq .  (12)  is  then  written, 

transformed  back  into  the  Schroedinger  representation,  as: 


0  -aw12p12  -awi3p13  -aw14p14 


[p]  = 


AW12P21 


°  -AW23P23  "AW24P24 


AW13P31  AW23P32  °  _AW34  P34 


AW14P41  AW24P42  AW34P43  ° 


2PU+P22+P33 

"3P12+P34 

-3P13“P24 

-2p14 

-3P21+P43 

-2P22+P11+P44 

-2P23 

~3P24_P13 

"3P31~P42 

_2P32 

-2P33+P11+P44 

-3P34+P12 

-2P41 

-3P42“P31 

_3p43+P21 

-2p44+p22+p 

where  the  first  matrix  in  Eq.  (13)  is  the  result  of  the  evaluation 
of  the  first  order  expansion  in  Eq.  (12),  with  =  (E^-Ej)/h, 
where  the  bar  represents  a  time  average  over  the  ij  collision. 

The  second  matrix  in  Eq.  (13)  is  taken  from  the  results  in 
Ref.  (5),  where  from  assumption  2  above,  we  have: 


1  /T  =  j  =  j  =  J  . 
Q  o  1  2 
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The  J^s  are  the  notation  in  Ref.  (5)  for  the  general  decay  rates 
due  to  a  quadrupolar  interaction.  The  condition  in  Eq.  (14) 
results  in  a  single  exponential  decay  for  the  nuclear  ensemble 
when  the  cell  is  oriented  at  the  magic  angle.  This  is  observed 
experimentally . 


The  observable  of  interest  in  this  particular  experiment  is 

1  12 

the  x-component  of  the  noble-gas  polarization  '  ,  <IX>  •  The 

corresponding  operator  for  the  spin  3/2  system,  written  in 
matrix  form  is: 


I 


x 


(1/2) 


0/300 
/3  0  2  0 

0  2  0  /3 

0  0/30 


The  expectation  value,  <Ix>,  is  given  by: 


<1x>  ■  Tr  [v  • 


(15) 


(16) 


For  the  spin  3/2  system,  Eq.  (16)  becomes: 


<Xx> 


/? 


(p12+p21) 


2(p23+p32) 


a  ( 
2  ' 


p43+p34) 


(17) 
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where  the  p^'s  are  the  elements  of  the  density  matrix,  with 


p^  the  population  density 
is  given  by  in  Eq.  (6) . 

of  the  state  whose 

With  the  following 

energy  eigenvalue 

definition: 

Mxl  = 

2  (p12  +  p21) 

(18a) 

Mx2  " 

(p23  +  P32) 

(18b) 

and , 

M  ,  = 

x3 

1  (p34  +  p43) 

(18c) 

we  can  write  Eq. 

(17)  as: 

<Jx>  " 

Mxl  +  Mx2  +  Mx3  * 

(19) 

We  have  adopted  the  symbol  'M'  for  the  sake  of  similarity  to  the 
2 

Bloch  equations  .  The  M's  here  represent  a  two-state  polariza¬ 
tion  of  the  four  state  spin  3/2  system.  In  exactly  the  same 
manner,  we  can  express  the  y-component  of  noble-gas  polarization 
in  terms  of  the  two  state  polarizations,  M^,  anc*  My3' 

appropriately  defined.  Using  Eqs.  (13  and  17-19)  it  is  only  a 
matter  of  algebra  to  show  that  the  following  set  of  equations 
describe  the  evolution  of  the  nuclear  ensemble  under  the 
quadrupolar  interaction  in  terms  of  the  two  state  polarizations: 


Xi 

-3/tq 

CM 

1^ 

< 

1 

0 

0 

1/tq 

0 

M  . 
xl 

V 

AW12 

-3/TQ 

0 

0 

0 

1/tq 

M  . 

yi 

**2 

0 

0 

-2/Tq 

-™23 

0 

0 

Mx2 

My2 

0 

0 

AW23 

-2/TQ 

0 

0 

My2 

A  .. 

x3 

1/tq 

0 

0 

0 

-3/tq 

-aw34 

Mx3 

A  , 

M 

0 

!/Tq 

0 

0 

4W34 

-3/Tq 

_My3_ 

(20) 
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Eq.  (20)  describes  the  evolution  of  a  set  of  three  oscillators. 
The  coupling,  which  results  from  the  inclusion  of  the  second 
order  terms,  is  between  the  first  and  third  oscillators.  If  we 
neglect  the  coupling  terms,  the  system  equations  are  easily  seen 
to  reduce  to  the  case  treated  in  Ref.  (1),  with  the  solution 
given  by  Eq.  (18)  of  Ref.  (1).  In  addition,  the  uncoupled 
equations  can  each  be  written  in  the  form: 


M  =  ( W  X  M)  +  M/Tq  (21) 

2 

which  is  exactly  the  form  of  the  Bloch  equations. 

Due  to  other  relaxation  effects  on  the  131Xe  ensemble,  we 

have  added  an  additional  phenomenological  term  to  all  the 

diagonal  elements.  This  procedure  agrees  with  the  results  we 

have  seen  for  the  situation  in  which  the  quadrupole  interaction 

1 3 

is  negligibly  small.  The  system  of  equations  for  the  two-state 
polarizations  then  become: 


(22) 
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where  1/T~  is  the  relaxation  rate  due  to  the  other  processes 

^  131 

present  in  the  cell,  for  instance,  spin  exchange  of  the  Xe 

nuclear  spin  with  the  rubidium  valence  electron. 

III.  EXPERIMENT 

The  experimental  apparatus  and  the  generation  of  the  signal 

from  the  noble  gas  nuclear  ensemble  has  been  previously  discussed 
14 

in  detail.  We  will  describe  here  some  of  the  salient  features 
of  the  experimental  set-up  for  the  sake  of  the  continuity  of 
this  report. 

The  experimental  apparatus  is  shown  schematically  in 

Figure  1.  The  heart  of  the  apparatus  is  the  experimental  cell, 

which  is  a  1-ml  Pyrex  sphere.  The  cell  is  sealed  on  a  vacuum 

system  with  an  excess  of  natural  rubidium  metal,  0.1  Torr 
129  131 

enriched  Xe,  0.4  Torr  enriched  Xe,  10  Torr  N_,  100  Torr  He 

129  ^  129 

and  10  Torr  H_ .  Enriched  Xe  consists  of  61.22  mole  %  of  Xe 

*  131  131 

and  13.44  mole  %  Xe.  Enriched  Xe  contains  63.89  mole  % 

131  129 

Xe  and  0.7  mole  %  Xe.  Remaining  components  of  each  of  the 

enriched  species  consist  of  even  isotopes  of  xenon. Since 

129  131  129 

there  is  a  Xe  impurity  in  enriched  Xe,  enriched  Xe  is 

129 

added  to  the  cell  to  ensure  a  sufficient  Xe  signal  so  that 

129 

the  effects  of  the  Xe  on  the  signal  can  be  properly  removed 
in  the  data  analysis.  The  ^  and  the  He  are  present  to  maximize 
the  rubidium  polarization,  which  is  used  to  polarize  the  Xe 
noble  gas  ensembles  through  spin  exchange.*4 

In  order  to  simplify  analysis  of  the  experimental  data  in 

the  light  of  the  theory  developed  in  Section  II,  the  electric 

131 

field  gradient  experienced  by  the  ensemble  of  Xe  nuclei  during 
collisions  with  the  cell  wall  must  be  nearly  cy lindrically 
symmetric.  It  was  found  in  the  course  of  this  work,  that  a  cell 
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I 

exhibiting  field  gradients  of  a  cylindrically  symmetric  nature 

could  be  made  by  adding  ^  to  the  cell  and  baking  the  cell  at 

80°C  for  about  two  weeks. 16  During  the  baking  period,  the 

rubidium  metal  was  held  in  the  small  region  of  the  tip-off  area 

of  the  spherical  cell  by  forcing  the  tip  region  to  be  cool 

relative  to  the  cell  walls.  Pyrex  walls  of  the  cells  so 

processed  were  covered  with  milky  particles,  which  we  later 

1 7 

identified  as  RbH.  Excess  rubidium  metal  was  still  present 

in  the  cell.  When  the  cell  wall  was  coated  uniformly  with  RbH, 

the  imaginary  axis  diametrically  through  the  tip  of  the  cell 

was  found  to  correspond  to  the  axis  of  cylindrical  symmetry  of 

the  quadrupolar  interaction.  This  condition  was  verified 

experimentally  if  a  cell  showed  the  magic  angle  behavior  when 

the  angle  between  the  magnetic  field  and  the  imaginary  axis 

through  the  tip  was  approximately  55°.  No  such  symmetry  could 

be  observed  in  the  cells  for  which  no  processing  was  done.  It 

131 

was  further  noted  that  Xe  polarization  is  much  higher  in  the 
RbH  coated  cells  than  in  the  cells  without  coating. 

Within  the  experimental  apparatus  the  cell  is  centered  in  a 
resistance-heated  oven,  provided  with  Pyrex  windows  at  either 
end.  The  oven  controls  the  temperature  to  within  ±0.1°C  with  a 
temperature  uniformity  across  the  cell  of  about  0.5°C.  The  cell 
is  mounted  in  a  holder  which  allows  it  to  be  rotated  in  two 
dimensions  from  outside  the  apparatus  with  a  precision  of  0.5 
deg.  The  oven  is  contained  in  a  cylindrical  coil  form,  providing 
three  mutually  perpendicular  Helmholtz  coils.  The  oven  assembly 
is  within  four  cylindrically  concentric  magnetic  shields.  These 
shields  reduce  the  external  magnetic  fields  to  below  10  uG. 

The  instantaneous  transverse  polarization  of  the  noble  gas 
ensembles  are  determined  by  measuring  the  effective  magnetic 
fields  of  these  polarized  ensembles.  The  nuclear  e»  'embles  are 
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polarized  in  spin  exchange  collisions  with  the  optically  oriented 
rubidium  vapor.  The  rubidium  vapor  is  optically  pumped  by 
illuminating  the  cell  with  circularly  polarized  light  from  a 
rubidium  discharge  lamp.  Detection  of  the  effective  magnetic 
fields  generated  by  the  polarized  nuclear  ensembles  is  accom¬ 
plished  by  utilizing  the  rubidium  vapor  as  a  magnetometer.  This 

18  19 

technique  has  been  thoroughly  discussed  in  the  literature  ' 

14 

and  our  particular  mechanization  has  been  previously  described. 
The  axes  for  this  work  have  been  labelled  in  a  slightly  different 
manner  than  what  was  done  in  Ref.  (14).  This  change  was  made  to 
facilitate  the  form  of  the  model  equations  and  presents  no 
difficulty  in  applying  the  results  as  displayed  in  Ref.  (14)  to 
this  work. 

The  particular  experimental  procedures  for  this  work  are  as 
follows.  In  the  presence  of  the  light,  a  2  mG  field  is  applied 
initially  along  the  y-axis,  as  defined  in  Figure  1.  This  condi¬ 
tion  continues  for  10  min.,  in  which  time  a  significant  noble 
gas  nuclear  polarization  is  built  up.  The  y-field  is  then 
switched  to  zero  while  a  precessional  magnetic  field  of  350  pG 
is  applied  along  the  z-axis  together  with  an  ac  magnetic  field 
(H  )  along  the  x-axis,  which  defines  the  sensitive  axis  of  the 
rubidium  magnetometer.  The  cell  is  arranged  in  the  apparatus 
with  the  cell's  imaginary  axis  in  the  x-z  plane.  The  imaginary 
axis  of  the  cell  can  be  rotated  a  full  360  degrees  in  this  plane. 
As  we  noted  earlier,  the  imaginary  axis  diametrically  through 
the  tip  of  the  cell  has  been  found  to  be  the  cylindrically 
symmetric  axis  of  the  quadrupole  interaction.  Rotation  of  the 
imaginary  axis  in  the  x-z  plane  then  corresponds  to  a  rotation 
of  the  cell  through  the  angle  as  defined  in  Eq.  (3).  There  is 
an  additional  ability  to  rotate  the  cell's  imaginary  axis  ±5 
degrees  out  of  the  x-z  plane.  This  feature  was  employed  to  check 
the  assumption  of  cylindrical  symmetry  of  the  interaction.  The 


2-16 


404577 


geometry  of  the  cell  is  shown  in  Figure  2.  4000  data  points  at 

equal  time  intervals  over  300  seconds  are  taken  after  the 
y-magnetic  field  is  switched  off.  Angle  data  at  10  deg.  cell 
angle  increments  in  the  x-z  plane  was  taken.  The  data  from 
the  first  half  cycle,  0  to  180  deg.  is  presented  in  this  report. 
All  data  in  this  report  were  acquired  at  a  cell  temperature  of 
8 1 . 3°C  ±0 . 5°C . 

IV.  ANALYSIS 

A.  The  signal  form.  In  order  to  facilitate  the  data  reduction, 
we  re-write  the  form  of  the  two-state  polarization  equations 
from  that  given  in  Section  II.  For  brevity,  we  introduce  the 
following  notation: 


w,  =  AW, .  =  w  3aa  (23d) 

3  34  c 

where  a  represents  a  time  average  of  a  over  the  collisions. 

Eq .  (22)  can  be  conveniently  broken  in  two  parts,  due  to  the 

particular  coupling  among  the  oscillators.  One  equation  to 
represent  the  second  oscillator,  subscripted  '2',  with  frequency 
oj 2  anc*  a  set  °f  equations  for  oscillators  '1'  and  '3',  with 
frequencies  and  respectively.  We  take  advantage  of  the 

symmetric  separation  in  these  two  oscillator  frequencies  from 
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the  center  frequency,  ui  ,  and  write  M  for  oscillator  1  and 

C  X 

3  as : 


M  =  M°  e-(P+r,t  F(8)  sin  (w  t) 
xs  xs  c 


F(e) 


cosh(et)  +  (r/c)  sinh(et)  r>3aa 
1  +  rt  r=3aa 
cos(et)  +  (r/e)  sin(et)  r<3aa 


2 

e 


I r 2  -  (3  aa)2| 


where  the  subscript  's'  denotes  either  1  or  3  and  Mxs  is  the 
initial  amplitude  of  the  signal  for  these  oscillators.  The 
frequency  of  the  second  oscillator,  can  be  written  as: 


where 


=  w  +  A  oi  _ 
2  c  2 


Aw. 


-22  -22 

108ct  D  +  21a  c 

w  4u 


(24) 


(25a) 


(25b) 


(25c) 


( 2  5d) 


(26) 


(27) 
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A u>2  is  due  to  a  second  order  effect,  which  is  negligibly  small. 
In  this  limit  we  can  write  the  expression  for  the  second 
oscillator  as: 


M 


x  2 


„o  -pt 
=  M  e  r 
xc 


sin  (wct) 


(28) 


If  we  assume  that  the  initial  amplitudes  of  the  first  and  third 
oscillators  are  equal,  we  can  then  write  the  expression  for  the 
signal  as  given  in  Eq.  (19)  as: 


Mx  =  Mjc  e  pt  [1  +  Me  rt  F(  0)  ]  sin  (^t) 


(29) 


where  the  ratio  of  the  oscillator  amplitudes,  M,  is  given  by: 


M 


2M°  /M 
xs 


o 

xc 


(30) 


At  the  magic  angle,  i.e.,  where  a  =  0,  Eq.  (29)  is  seen  to  reduce 
to : 


M  =  M°  [1  +  M]  e  sin  (w  t) 

X  xc  c 


(31) 


Eq.  (31)  is  seen  to  represent  a  single  exponentially  decaying 
sinusoid,  which  is  what  is  experimentally  observed.'*' 

131 

B.  Data  reduction.  As  we  previously  stated,  Xe  is  not 

129 

commercially  available  without  a  Xe  impurity.  For  this 

12  9 

reason,  we  chose  to  enhance  the  unwanted  Xe  signal  by  adding 
129 

enriched  Xe  to  the  cell  so  that  we  could  ensure  its  complete 
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removal  from  the  data  before  we  analyzed  the  decay  of  the  Xe 
ensemble . 

Typical  raw  data  is  displayed  in  Figure  3.  This  data  was 

taken  off  the  magic  angle,  thus  the  non-exponential  character  of 

the  ^*Xe  decay  is  clearly  observed.  Since  the  nuclear  spin  of 
129 

Xe  is  1/2,  its  decay  at  all  cell  orientations  is  simple 

12  9 

exponential.  The  form  of  the  Xe  decay  is  thus  exactly  known 

and  we  then  utilized  a  Kalman  filtering  technique  to  remove  this 

signal  from  the  raw  data.  The  Kalman  technique  was  used  since 

this  method  was  designed  to  remove  'signals'  buried  in  noise. 

131 

In  this  situation,  we  simply  treat  the  Xe  signal  and  the 

noise  as  'noise'  and  remove  a  single  exponentially  decaying 

129 

sinusoid  representing  the  effects  of  Xe.  The  result  of  this 

processing  is  displayed  in  Figure  4.  A  'fast  fourier  transform' 

129 

was  then  performed  on  the  data  with  the  Xe  signal  removed  to 
129 

verify  that  the  Xe  had  been  completely  removed. 

Due  to  the  subtle  differences  between  data  fitting  with  a 

22 

Kalman  technique  and  a  Least  Squares  method  ,  we  chose  to  employ 

131 

a  Least  Squares  routine  to  fit  the  Xe  data  to  our  model.  We 
used  an  eight  parameter  fitting  routine  of  the  form: 

-(x.+2xc-)t  -x_  t 

S  =  x3e  [1  +  Me  F(e)]  sin  (x1+x2t)  +x7+xgt  (32) 

cos(x,t)  +  (xc/x,)  sin  (x,t)  (33a) 

F (e)  =  656  6 

cosh  (Xgt)  +  (x^/Xg)  sinh  (x^t)  ,  (33b) 

Eq.  (32)  was  used  with  Eq.'s  (33a)  and  (33b)  for  the  angle  data 
far  from  the  magic  angle  and  very  near  the  magic  angle. 
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respectively . 

represent  the 

In  Eq.'s  (32  through  33b), 

following  terms: 

the  eight  parameters 

X1  = 

initial  phases 

(34a) 

X2  = 

uc 

(34b) 

X3  = 

M° 

xc 

(34c) 

Xj  =  dc  offset  of  the  signal 


x„  =  linear  drift  in  the  dc  offset. 


( 34d) 

(34e) 

(34f) 

(34g) 

( 34h) 


The  factor  t./(t  +t„)  in  Eq.  ( 34 f )  results  from  the  time 

S  S  v 

average  of  the  dephasing  terms  over  the  duration  of  the  collision. 
ts  is  the  duration  of  the  interaction,  that  is  the  sticking  time 
on  the  wall,  and  xv  is  the  time  between  collision. 

In  the  fitting  routine,  we  have  set  the  initial  phases  of 
the  three  oscillators  to  be  equal.  There  is  no  a  priori  reason 
why  this  should  not  be  and  we  observe  nothing  in  the  results  of 
our  data  fit  to  indicate  that  this  is  not  a  valid  assumption. 
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We  further  had  to  make  an  assumption  concerning  the  ratio  of  the 

oscillator  amplitudes,  M  given  by  Eq.  (30).  An  attempt  to  let 

this  vary  in  the  Least  Squares  routine  resulted  in  instabilities 

in  the  data  fit.  In  Ref.  (1) ,  we  considered  the  case  of  low 
2  3 

polarization.  In  this  limit,  using  Eq.  (8  and  10)  of  Ref.  (1) 
and  Eqs.  (18a-c)  we  find: 


3:4:3 


(35) 


In  the  limit  of  complete  polarization,  on  the  other  hand,  which 
implies  p^^l  in  Eq.  (8)  of  Ref.  (1),  we  find: 


1:2:1 


(36) 


We  have  good  reason  to  believe  that  for  these  experiments,  the 
nuclear  ensemble  was  not  in  the  low  polarization  limit,  however, 
we  also  definitely  know  that  the  nuclear  ensemble  did  not  attain 
100%  polarization,  this  implied: 

1  <  M  <  1 . 5  (37) 

where  1  is  the  limit  of  complete  polarization  and  1.5  is  the  low 
polarization  limit.  We  additionally  demand  that  this  ratio 
remain  constant  over  all  the  angles  at  a  given  temperature  and 
pumping  time.  This  is,  as  was  argued  in  Ref.  (1),  due  to  the 
fact  that  the  first  order  quadrupolar  interaction  does  not  enter 
into  the  longitudinal  relaxation  of  the  noble  gas  ensemble  and 
hence  cannot  affect  the  initial  amplitude  of  the  signal.  By 
employing  different  ratios  under  the  above  constraints  in  the 
eight  parameter  fitting  routine,  we  found  that  M  =  1.35  gave  the 
best  results.  This  value  was  then  fixed  in  the  routine. 
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V.  RESULTS 


The  three  parameters  of  interest  in  the  eight  parameter 
Least  Squares  routine  are  x^  ,  and  x^ .  From  these  terms  we 
derive  the  phenomenological  relaxation  rate, 

~  =  (3.76  ±  .41) x  10-3  sec-1  (38) 

2 

and  the  transverse  quadrupole  relaxation  rate, 

=  (5.94  ±  .65)  x  10~4  sec-1  .  (39) 

TQ 

The  angularly  dependent  dephasing  rate,  3¥a ,  is  displayed  in 
Figure  5.  The  errors  associated  with  the  derived  terms  represent 
one  standard  deviation.  The  solid  line  in  Figure  5  represents 
a  Least  Squares  fit  of  the  time  average  of  the  dephasing  rate  to 
the  equation: 


A  ( 0 )  =  j  [3  cos2  6  -  1]  .  (40) 

where  9  is  the  angle  between  the  magnetic  field  axis  (z-axis) 
and  the  imaginary  axis  diametrically  through  the  tip  of  the  cell. 
From  the  fit  we  find: 

A ( 0 )  =  (8.9  ±  0.4)  v  10~3  [3  cos2  9  -  1)  .  (41) 

Comparing  Eq .  (40)  and  Eq .  (41)  we  find: 

j  e2Qq  =  3.56  x  10-2  sec-1  .  (42) 


^Ts  +  T 
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We  can  estimate  the  sticking  time  on  the  wall, 
simple  adsorption  model.  The  sticking  time  is 


t  ,  by  assuming  a 

s  24 

then  given  by  : 


=  Toe 


E  /kT 

a 


(43) 


energy 
i  «-ll 

10  sec. 


-12 

sec,  and  E  is  the  adsorption 
13  a 

Taking  E  to  be  0.13  eV  ,  we  compute:  t  =  6.6  x 

d  S 


where  t  is  taken  to  be  10 
o 


t  is  the  time  between  wall  collisions,  and  it  can  be 

V  25 

estimated  from  simple  gas  kinetic  arguments: 


t  =4  Vol/vA  (44) 

V  s 

where  Vol  is  the  volume  of  the  cell,  Ag  is  the  surface  area  and 
v  is  the  kinetic  velocity  of  the  xenon  atoms.  Evaluating 

Eq.  (44)  for  our  experimental  conditions  yields:  t  =  3.5  x 

-  5  v 

10  sec.  Since  t  < < t  ,  we  can  neglect  t  in  the  denominator 

S  V  s 

of  Eq.  (42)  and  calculate  the  average  electric  field  gradient: 

eq  =  -3.45  x  1G^3  statcoulombs  cm  3  .  (45) 

26 

Recently  Keenan  et  al  estimated  the  electric  field  gradient 

131  . 

which  arises  at  the  nucleus  of  a  Xe  atom  due  to  a  distortion 
of  the  spherical  shells.  This  distortion  takes  the  form  of 
electronic  excitation  with  dipole  selection  rules.  They  assumed 
a  one  electron  model  and  the  excitation  5p^  -*  5p56s.  The  field 
gradient  was  then  estimated  to  be: 

i  o  —3 

eq  =  -7.0  x  10  statcoulombs  cm  .  (46) 

Given  the  approximations,  we  find  our  value  in  good  agreement 
with  this  estimate. 
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The  transverse  quadrupole  rate  is  given  in  Ref.  (5)  as: 


36n 

5 


TS 


T  +T 
S  V 


e2Qq 

3 


T 


C 


(47) 


From  our  measured  value  of  the  quadrupole  relaxation  rate  given 

in  Eq .  (39),  and  our  estimate  of  the  electric  field  gradient 

strength  given  in  Eq .  (45) ,  we  compute  the  correlation  time  for 

this  interaction  to  be:  t  =  3.5  x  10-7  sec. 

c 


Finally  as  a  check  on  our  parameter  fit,  we  compare  the 

131 

phenomenological  relaxation  rate  of  the  Xe  with  the  decay 

rate  of  the  129Xe  signal.  From  previous  measurements13,14  on 
129  131 

the  decay  of  Xe  and  Xe,  we  predict  that  at  this  temperature 
the  predominant  relaxation  is  due  to  spin  exchange  collision  of 
the  xenon  atoms  with  the  rubidium  vapor.  We  then  write: 


NRb  CTex  Vrel 


(48) 


where  is  the  rubidium  density,  is  the  Rb-Xe  spin 

exchange  cross  section  and  V  ^  is  the  relative  Rb-Xe  velocity. 

The  ratio  of  relaxation  rates  is  simply  the  ratio  of  cross 

.  129 

section  for  spin  exchange.  From  the  Xe  signal,  which  we 

extracted  from  the  original  data,  we  found: 

1/T  =  (l .  34  ±  0.7)  x  10-2  sec-1  (49) 

and  thus  we  estimate  the  cross  section  ratio  to  be: 

129  131  .  /tec 

0  /a  1/3.56  .  (50) 

ex  ex 
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From  Herman's  model  of  Rb-noble-gas  spin  exchange  we  predict  the 

2  7 

ratio  of  cross  sections  to  be  : 

129  111 

a  / a  =  1/2.28  .  (51) 

ex  ex  '  ' 

We  believe  that  this  is  very  good  agreement  supporting  the 
validity  of  our  fitting  procedure,  and  that  the  difference 
between  Herman's  prediction  and  our  estimate  is  due  to  the 
contribution  of  the  Rb-Xe  van  der  Waals  molecules  in  the 
alkali-noble-gas  spin  exchange  interaction. ^ 


VI.  DISCUSSION 

We  have  investigated  the  nuclear  transverse  spin  relaxation 

of  13^Xe  under  a  quadrupole  coupling  between  the  *  Xe  nucleus 

and  the  electric  field  gradients  which  arise  at  the  nucleus 

during  collisions  of  the  xenon  atoms  with  the  walls  of  the  cell. 

We  picture  the  collisions  as  resulting  in  a  deformation  of  the 

xenon  electronic  shells.  It  is  this  deformation  which  produces 

the  electric  fields  at  the  xenon  nucleus.  The  angular  dependence 

of  this  interaction  is  interpreted  to  mean  that  adsorption  sites 

on  the  cell  walls  have  a  non-isotropic  distribution.  Further 

supporting  this  view  is  our  finding  of  a  dependence  of  this 

interaction  on  the  distribution  of  the  metalic  rubidium  on  the 

cell  walls.  A  nearly  uniform  distribution  of  the  rubidium 

reservoir  over  the  surface  of  the  cell  results  in  a  very  strong 

131 

relaxation  of  the  Xe  ensemble,  in  fact  complete  relaxation  is 
found  to  take  place  in  a  few  precessional  periods,  and  thus 
usable  data  cannot  be  obtained.  In  the  cells  that  are  uniformly 
coated  with  RbH,  the  adsorption  sites  on  the  walls  appear  to  be 
isotropic  with  the  location  of  excess  rubidium  metal  being 
different  from  the  rest.  Relaxation  data  that  fit  to  our  model 
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are  obtained  only  in  those  RbH-coated  cells  in  which  the  rubidium 
metal  is  driven  to  the  tip  of  the  cell. 

We  believe  that  our  technique  of  observing  this  effect  is 
quite  valuable  since  we  obtain  information  from  both  the 
dephasing  and  transverse  relaxation,  which  have  different 
dependencies  on  the  interaction  parameters.  We  have,  however, 
still  not  been  able  to  definitely  determine  the  adsorption 
energy  or  sticking  time  of  the  xenon  atoms  on  the  cell  wall  or 
the  correlation  time  of  the  interaction,  and  thus  a  precise  value 
of  the  electric  field  gradient,  which  we  estimated,  remains  in 
some  question.  We  have  begun  experiments  to  investigate  these 
two  parameters  which  will  allow  us  to  make  a  better  determination 
of  the  interaction  strength. 
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FIGURE  CAPTIONS 

Figure  1.  Schematic  representation  of  the  experimental 
apparatus 

Figure  2.  Geometrical  arrangements  of  the  experimental  cell 

Figure  3.  Typical  example  of  the  raw  signal.  Low  frequency 

131 

component  is  due  to  the  Xe  precession,  while  the 

129 

high  frequency  is  the  Xe  signal.  This  data  was 
taken  at  the  cell  orientation  at  which  the  magnetic 
field  axis  coincides  with  the  imaginary  axis 
diametrically  through  the  tip  of  the  spherical  cell. 

Figure  4.  Typical  example  of  the  signal  after  the  effects  of 
129 

the  Xe  have  been  removed.  The  data  shown  in 
Figure  2-3  are  displayed. 

Figure  5.  Quadrupole  dephasing  terms  as  a  function  of  angle. 

The  error  bars  on  the  data  points  represent  one 
standard  deviation  and  the  solid  line  is  a  Least 
Squares  fit  to  the  data. 
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SECTION  III 

NOBLE  GAS  NUCLEAR  FREQUENCY  SHIFT 

3 . 1  INTRODUCTION 

The  shift  of  the  precessionai  frequency  of  the  noble  gas  nucleus 
in  a  magnetic  field  reported  in  this  section  is  associated  with 
the  transfer  of  angular  momentum  from  an  electronically  polarized 
alkali  atom  to  an  unpolarized  noble  gas  nucleus  in  the  spin 
exchange  process.  Such  frequency  shift  has  been  reported  for  the 
case  of  free  electrons  interacting  with  optically  pumped  alkali 

atoms1  1 .  No  such  information  is  yet  available  to  our  knowledge 

83  129  131 

for  noble  ga.'*  nuclei,  such  as  Kr  ,  Xe  and  Xe  . 

Since  the  angular  rate  in  the  nuclear  magnetic  resonance  (NMR) 
gyro  is  sensed  as  a  shift  in  the  precession  frequency  of  certain 
nuclear  species,  the  information  regarding  the  frequency  shift 
can  lead  to  a  determination  of  gyro  bias  shift  due  to  the  spin 
exchange  interaction  with  the  polarized  alkali-vapor.  As  will  be 
shown  later,  the  frequency  shift  of  the  noble  gas  nucleus  is 
strongly  dependent  on  both  the  number  density  of  the  alkali-atom 
and  the  alkali  polarization.  This  would  mean  that  the  same 
information  would  also  lead  to  determinations  of  temperature 
dependence  and  the  alkali-pumping  light  intensity  dependence  of 
the  NMR  gyro  bias. 

In  this  section  we  review  briefly  the  theoretical  background  of 
the  spin  exchange  interaction  to  give  an  understanding  of  the 
frequency  shift,  and  subsequently  present  the  experimental 
results. 
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3.2  THEORY  OF  FREQUENCY  SHIFT 

In  order  to  review  the  theory,  we  will  rely  on  the  results 

4 

previously  reported  ,  and  repeat  only  selected  details  here  as 
clarity  demands. 

We  start  with  Herman's  effective  interaction  Hamiltonian  between 

5 

an  alkali  atom  and  noble  gas  nucleus  in  the  form  : 

3Ceff  =  (pp)  (52) 


where  I  is  the  noble-gas  nuclear  spin,  S  the  alkali  electronic 
spin,  gN  the  nuclear  g-f actor,  the  nuclear  magneton,  and 

the  Bohr  magneton.  V  is  the  exchange  enhancement  factor,  given 
by . 


^  (R) 


i  .  £  <t|i>  ui(B) 

i=2  UjTB) 


(53) 


where  R  is  the  internuclear  separation  between  the  alkali  and 

noble  gas  atoms.  The  summation  over  i  refers  to  the  alkali 

valence  electron  (1)  and  the  noble  gas  electrons  (2,3...,n). 

The  factors  <i|l>  are  the  overlap  integrals  and  u^(B)  is  the 

amplitude  of  the  ith  orbital  of  the  isolated  atom  at  the  noble 

gas  nucleus,  B.  For  the  internuclear  separation  of  R  =  7.8a 

“8  ° 

(a  is  the  unit  of  atomic  distance;  a_  =  0.529X10  cm),  the  values 
o  °  4 

of  v  for  the  Rb-Kr  and  Rb-Xe  system  have  been  calculated  to  be  : 

r?Kr  "  40  and  ^xe  =  “59  •  (54) 
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We  consider  a  gaseous  mixture  of  alkali  and  noble  gas  atoms 
subjected  from  the  outside  only  to  a  constant  uniform  magnetic 
field,  but  interacting  with  each  other  through  the  relatively 
short  range  forces,  so  that  the  otherwise  free  atoms  occasionally 
collide  with  each  other.  These  collisions  may  be  either  of  the 
binary  type  or  molecular  forming  (three  body)  type.  All  colli¬ 
sions  are  such  that  the  mean  time  between  collisions,  r  ,  is  long 
compared  to  the  duration  of  the  collision,  td.  It  is  further 
assumed  that  the  external  magnetic  field  is  so  weak  that  the  pre¬ 
cession  frequencies  of  the  electron  moment  and  the  nuclear  moment 
are  orders  of  magnitude  smaller  than  the  strength  of  the  inter¬ 
action  given  by  Eq.  (52),  when  the  internuclear  separation  is 
within  the  atomic  range. 

In  what  follows  we  take  the  total  Hamiltonian  acting  on  the  spin 
variables  to  be: 


3C  =  -hq _S^  -  tioi  I  +  [h  y  (R)  7- s’]  (55) 

o  z  I  z 

where  the  first  two  terms  of  the  right  hand  side  of  the  equation 
are  the  coupling  to  the  external  magnetic  field,  which  will  be 
represented  by  K  ,  and  the  last  term  represents  the  exchange 
interaction,  where  y(R)  is  the  strength  of  the  interaction: 

V(R)  =  (i|i)  gN^Qu1(B)2  n{ R)2/h  .  (56) 


v(R)  falls  off  rapidly  with  internuclear  separation,  R,  and  thus 
the  kinematics  of  the  interaction  is  to  picture  the  spins  before 
a  collision  as  acted  on  by  JCQ,  then  coming  under  the  combined 
action  of  3CQ  +  hy(R)7* 7 during  the  interaction  and  finally,  as  the 
particles  separate,  going  back  to  the  uniform  precession  of  J(Q. 
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We  now  consider  a  much  simplified  system  consisting  of  an  atom 
with  a  single  electron  in  a  state  and  no  nuclear  spin  (the 

alkali  atom  S= 1/2)  and  an  atom  with  no  electronic  spin,  S=0  and 
nuclear  spin  1=1/2  (the  noble  gas  atom).  The  system  is  described 
quantum  mechanically  by  a  density  matrix  of  (2S  +  1)(2I  +  1) 
dimensions,  namely,  by  a  4  x  4  matrix. 

The  density  matrix  can  be  constructed  in  the  representation  of 

2  2 

the  ortho-normal  set  of  eigenfunctions  of  I  ,  Iz,  and  S  ,  Sz  for 
the  free  (noninteracting)  atoms.  The  time  development  of  the 
density  matrix  resulting  solely  from  the  spin  exchange  interac¬ 
tion  can  then  be  investigated  in  the  interaction  representation. 


The  derivation  of  the  shift  of  the  noble  gas  nuclear  precession 
frequency  involves  the  equation  of  motion  for  the  mean  value  of 

the  operator  I+  (or  I  ),  where,  as  commonly  defined,  I*  =  I  ±  i  Iy, 

••  X  6 
in  the  Schrodinger  representation.  Following  the  work  of  Major  , 

it  has  been  shown  that 


<I+>  = 


,  ["(i/2)y2  i-j? 

1  ave  D 

ta  l+y2  rj? 

L  ave  D 


-  i  <S  > 
z 


yaveTD 


<I+> 


> 


where  7  is  the  strength  of  the  spin  exchange  interaction, 
3VG 

given  by  Eq.  (56),  averaged  over  the  collisional  parameters. 
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The  real  part  of  the  coefficient  of  <I+>  represents  the  phase 
coherence  decay  time,  Tex,  due  to  the  spin-exchange  interaction, 
i.  e.  , 


-1 


ex 


■  ^ 


^aveTD 


y2  Tn 
ave  D 


(58) 


The  imaginary  term,  which  is  seen  to  be  proportional  to  the  elec¬ 
tron  spin  polarization,  represents  a  shift  in  the  precession  fre¬ 
quency  of  the  nuclear  spin  in  the  external  magnetic  field: 


A  Cd 


y 

ave  D 


1  + 


y2  t2 
ave  D 


<S  > 

z 


(5Q) 


This  shift  can  be  written  as  an  equivalent  magnetic  field,  given 
by: 


AH 


y  T 

ave  D 


1  + 


y2  t2 
ave  D 


<S  > 

z 


(60) 


where  AH  is  now  a  magnetic  field  'experienced'  by  the  noble  gas 
atoms  due  to  the  presence  of  the  alkali  atom. 


3.2.1  Frequency  Shift  Due  to  Exchange  Process  of  the  Binary  Type 

Consider  the  spin  exchange  process  due  to  the  binary  type  colli¬ 
sions,  for  which  V  «  l7.  Eq.  (60)  can  now  be  written  as: 

'  ave  D 


A  H 


y 

ave 


<S  > 

z 


(61) 
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The  collision  rate  can  be  written  as: 


7 D  -  ’ d 2  NA 


rel 


(62) 


where  d  is  the  mean  interaction  distance  and  N»  is  the  alkali 
atom  number  density.  vre^  is  mean  relative  velocity  between 

the  alkali  atoms  and  noble-gas  atoms.  The  binary  collision  dura¬ 
tion  is  taken  to  be: 


T 

D 


d 

vrel 


(63) 


Combining  Eq's.  (62)  and  (63)  with  the  definition  of  y  given  in 
Eq.  (56),  Eq.  (61)  is  found  to  be: 


AH 


Na^o  u1(B)2  4(R)2 


<S  > 


z 


(64) 


where  we  have  appropriately  averaged  values  of  u^(B)  and  4(R). 
Since  u.(B)  is  the  wave  function  of  the  unperturbed  alkali 
valence  electronic  orbital,  2jrdJ  u^  (B)  is  of  order  unity. 

We  note  in  Eq.  (64)  that  the  value  of  AH  for  different  isotopes 
of  the  same  noble  gas  atom  is  the  same  for  the  same  alkali  polar¬ 
ization  environment  (neglecting  a  small  effect  of  the  different 
nuclear  structure),  but  that  different  noble  gas  species  experi¬ 
ence  different  magnetic  field  shifts,  mainly  due  to  the  different 
values  of  the  exchange  enchancement  factor,  tj  . 
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If  we  pick  <Sz>  -  0.1  and  a  sample  temperature  of  65°C,  we  find 
the  magnetic  field  shifts  for  Kr  and  Xe  due  to  the  presence  of 
Rb  as : 


AH  Kr  ~  2.6  ^G,  and 
AH  Xe~  5,7  nG  . 


(65) 


Under  the  present  mechanization  of  Litton 's  NMR  gyro,  the  gyro 

bias  shift,  Q  ,  due  to  the  magnetic  fields  seen  by  the  noble  gas 
a  8 

atoms  is  shown  to  be  : 


a 


(66) 


where  the  subscripts  a  and  b  refer  to  the  two  nuclear  species  to 

be  used  in  the  sample,  and  the  Y's  are  the  corresponding  nuclear 

gyro  magnetic  ratios.  When  the  results  of  Eq.  (65)  are  used,  the 

gyro  bias  shift  becomes  580  deg/hr.  Samples  that  contain  a  pair 

129 

of  isotopic  nuclei  of  the  same  noble  gas  atom,  such  as  Xe  and 
Xe13\  will  show  essentially  zero  bias  shift  since  the  factor 
(.AH^  -  AH^ )  is  diminished. 

In  an  NMR  gyro,  this  gyro  bias  shift  should  not  be  a  problem  as 
long  as  it  remained  constant.  Eq.  (64)  indicates,  however,  the 
gyro  bias  shift  will  vary  as  a  function  of  the  cell  temperature 
(temperature  sensitivity)  or  of  the  alkali-pumping  light  inten¬ 
sity.  These  shifts  can  be  clearly  deduced  by  noting  that  a 
change  of  the  cell  temperature  will  result  in  a  change  of  both 
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the  density  of  alkali  atoms,  N.  ,  and  the  alkali-polarization,  S 

and  that  a  variation  of  the  pumping  light  intensity  at  a  fixed 

temperature  will  result  in  a  variation  of  <  S  -> ,  It  is  for  this 

z 

reason  that  two  nuclear  isotopes  of  the  same  noble  gas  atoms, 

129  131 

i.e.  Xe  and  Xe  ,  are  the  preferred  embodiment  in  the  Litton 
approach  to  the  development  of  an  NMR  gyro. 

3.2.2  Frequency  Shift  in  a  Pair  of  Isotopic  Nuclei  of  the  Same 
Noble  Gas 

We  have  argued  in  the  previous  section  that,  when  the  spin 
exchange  is  dominated  by  binary  type  collisions,  a  pair  of  iso¬ 
topic  nuclei  of  the  same  noble  gas  atom  experiences  the  same  fre¬ 
quency  shifts,  which  results  in  a  "zero"  gyro  bias  shift.  Inter¬ 
actions  other  than  the  dipolar  spin-exchange  process  that  may 
lead  to  a  frequency  shift  are  discussed  in  the  following. 

Nuclear  Structure  Effect 

There  will  be  a  limit  to  the  condition  that  a  pair  of  isotopic 
nuclei  of  the  same  noble  gas  atom  experiences  the  same  frequency 
shift.  This  is  due  to  the  isotopic  "hyperfine  anomaly"  resulting 
from  nuclear  structure  effects,  which  makes  the  strength  of  the 
hyperfine  interaction,  1(R) ,  not  strictly  the  same  for  each  mem¬ 
ber  of  the  pair.  Nevertheless,  the  departure  due  to  this  effect 
of  the  small  variance  in  the  gyromagnetic  ratios  of  the  pair  is 
never  more  than  one  part  in  10^. 

Effect  of  Molecular  Formation 

7  9 

It  has  been  shown  previously  '  that,  in  addition  to  the  binary 
collisions,  the  spin  exchange  takes  place  through  the  formation 
of  weakly-bound  molecules,  such  as  Rb-Xe,  in  the  presence  of  a 
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third  colliding  body,  such  as  ^  buffer  gas  (three  body  collisions 
If  the  duration  of  the  interaction  is  taken  to  be  the  duration  of 
the  three  body  collision,  r  ,  the  magnetic  field  experienced  by 
the  noble  gas  nuclei  due  to  the  formation  of  molecules,  AH 
shown  to  be^ : 


mole' 


is 


AH 


mole 


gNPN 


tf) 


ave 


1+1 


2 

ave 


<S  > 

z 


(67) 


where  T0  is  the  molecular  formation  rate,  and  both  T_  and  r  are 
r  r  c 

functions  of  the  buffer  gas  density.  It  remains  to  be  investi¬ 
gated  if  the  interaction  strength  and  the  collision  time  are  the 
same  for  different  isotopes  of  the  same  noble  gas  under  conditions 
of  identical  buffer  gas  density. 

Effect  of  Nuclear  Quadrupole  Moment 

131  — 

Consider  the  nucleus  of  Xe  (I  =  3/2),  the  nuclear  quadrupole 
moment  of  which  interacts  with  the  electric  field  gradient  at  the 
nuclear  site.  It  has  been  stated  in  Section  II  that  the  electric 
field  gradient  arises  from  deformation  of  the  otherwise  spherical 
electronic  shells  of  the  zenon  atom  during  inelastic  collisions 
with  the  cell  walls.  (The  quadrupole  interaction  is  absent  in  the 
nucleus  of  Xe  ,  because  I  =  1/2.)  The  quadrupole  interaction 
of  the  noble  gas  nucleus  is  modeled  as  21  oscillators,  namely 
3  oscillators  for  Xe^J^,  each  with  a  single  frequency.  The 
frequencies  are  dependent  on  the  orientation  of  the  cell  with 
respect  to  the  magnetic  field.  When  a  cell,  where  the  nuclear 
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quadrupole  interaction  is  made  cylindrically  symmetric,  is 
oriented  to  the  'magic'  angle,  the  three  oscillator  frequencies 
become  nearly  identical  to  the  frequency,  u>c,  given  by  Eq.  (23a)- 


OJ  = 

w  + 54ha 

(68a) 

c 

U) 

where 

U)  - 

-YhH 

(68b) 

b  = 

sinG  cose  . 

(68c) 

where  w  is  the  unperturbed  nuclear  precession  frequency  in  a 
magnetic  field  H,  and  0  is  the  angle  between  the  magnetic  field 
axis  and  the  axis  of  cylindrical  symmetry.  The  strength  of  the 
quadrupole  interaction,  a,  has  been  shown  to  be: 


a 


e2Qq 

41(21-1) 


(69) 


where  eq  and  eQ  are  the  electric  field  gradient  and  the  nuclear 

quadruDole  moment,  respectively,  as  described  in  Section  II. 

is  the  sticking  time  of  xenon  atoms  on  the  cell  walls,  and 

r  is  the  time  between  wall  collisions, 
v 

Any  departure  of  the  nuclear  precessional  frequency,  u>  ,  from  the 
unperturbed  frequency,  w,  results  in  the  frequency  shift.  The 
shift  is  seen  from  Eq.  (68a)  to  be  proportional  to  the  square 
of  the  strength  of  the  nuclear  quadrupole  interaction.  In  a 
high  magnetic  field  (w  is  large),  however,  this  frequency  shift 
could  be  made  relatively  small.  Discussions  so  far  have  been 
limited  to  the  first  order  effect  of  the  nuclear  quadrupole  inter¬ 
action.  The  second  order  effects  are  small  and  are  neglected. 
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3 . 3  EXPERIMENT 

3.3.1  Description  of  Experiment 

Shown  in  figure  6  is  the  experimental  arrangement  designed  to 
measure  the  magnetic  field  experienced  by  the  noble  gas  nucleus 
due  to  the  spin  exchange  interaction  with  an  optically  oriented 

7 

alkali  vapor.  This  setup  has  been  described  in  detail  elsewhere 
Modifications  of  the  arrangement  are  made  so  that  the  alkali 
pumping  light  beam  can  be  made  left  or  right  hand  circularly 
polarized  by  placing  suitable  polarizers  in  front  of  the  light 
source.  The  left  and  right  hand  circular  polarizers  are  made  by 
laminating  a  sheet  of  plastic  linear  polarizer  and  two  sheets  of 
plastic  retarder.  The  two  polarizers  so  made  transmit  nearly 
equal  intensity  of  rubidium  D1  light.  The  amount  of  circularly 
polarized  light  was  measured  to  be  92%  for  one  and  96%  for  the 
other  of  the  transmitted  light. 

The  experiments  were  performed  by,  first,  polarizing  the  noble 
gas  nuclear  spins  through  spin  exchange  collisions  with  the 
optically  oriented  rubidium  vapor  (Pump  Mode).  The  rubidium 
atoms  were  oriented  by  the  absorption  of  left  or  right  hand  cir¬ 
cularly  polarized  rubidium  light  in  a  standard  longitudinal 
pumping  scheme:  This  configuration  is  shown  in  figure  7(a). 

The  pumping  takes  place  in  the  presence  of  a  longitudinal  mag¬ 
netic  field  (Pump  Field)  of  approximately  100  mG  along  the  z-axis 
The  longitudinal  noble  gas  polarization,  I ,  has  been  shown  to 
obey  the  usual  pumping  law,  given  the  initial  condition 

■  i  ( t=0 )  --  =  <  I  °>: 
z  z 


<Iz(t)>  = 


<1 


|  .l-exp(-t/Tp)  J 


(70) 
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Figure  6.  Schematic  representation  of  the 
experimental  apparatus 
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where  is  the  "pumping"  time  of  the  ensemble  of  noble  gas 
nuclear  spins.  The  ensemble  was  pumped  for  a  time  sufficient  to 
build  up  a  significant  polarization.  This  time  can  range  from 
5  minutes  at  low  sample  temperatures  to  30  minutes  at  high  sample 
temperatures . 

When  a  sufficient  nuclear  polarization  was  built  up  along  the 
light  beam  axis  (z-axis),  the  pump  field  was  turned  off,  and,  at 
the  same  time,  a  "classic"  2  pulse  was  applied  along  the  x-axis 
to  the  nuclear  system.  This  "steers"  the  nuclear  polarization 
away  from  the  z-axis  into  the  x-y  plane  as  shown  in  figure  7(b). 

Finally,  the  precession  of  the  noble  gas  nuclear  spin  polariza¬ 
tion  about  the  z-axis  was  accomplished  by  applying  a  precession 
field  (~100  nG)  along  the  z-axis,  figure  7(c).  The  optically 
oriented  rubidium  vapor  in  this  scheme,  then,  acts  either  to 
directly  increase  or  decrease  the  effective  precession  field  as 
experienced  by  the  noble  gas  ensemble,  depending  on  the  circular 
polarizers  used  in  the  measurement. 


The  precession  frequency  of  the  ensemble  of  the  noble  gas  nuclei 
is,  correspondingly,  increased  or  decreased  such  that: 


w1  =  -V(Hq  +  AH) 

w-  =  -Y(H  -  AH) 
2  o 


(71) 


where  Hq  is  the  magnetic  field  applied  externally  along  the  pre¬ 
cession  axis  and  AH  is  the  field  experienced  by  the  noble  gas 
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ensemble  due  to  the  spin  exchange  interaction.  V  is  the  nuclear 

gyromagnetic  ratio  of  the  noble  gas  atoms.  The  subscripts,  1  and 

2,  refer  to  the  two  different  circular  polarizers.  From  these, 

AH  is  readily  obtained  by  measuring  the  precession  frequencies, 

and  u>  .  A  similar  experimental  arrangement  has  been  employed 
1  1  1-3 

previously  by  Balling,  et  al  ,  for  the  measurement  of  fre¬ 
quency  shifts  of  electrons  interacting  with  the  polarized  alkali 
vapor  through  the  spin  exchange  process. 

In  order  to  observe  the  precessing  nuclear  spin  polarization  sig¬ 
nal,  the  rubidium  vapor  was  mechanized  as  a  magnetometer.  The 
magnetometer  detects  magnetic  field  variations  along  the  x-axis 
due  to  the  precessing  noble  gas  nuclear  ensemble.  The  form  of 
the  observed  signal  has  been  shown  to  be7: 

Si  =  Saexp  f-t/(T2)alsin(u;ai  +  °ai) 

(72) 

tS^exp  t/(T2  )faj  sin(u>b^  +  0^ ) ,  for  i=l  and  2 

where  the  subscripts  a  and  b  refer  to  each  species  of  noble  gas 
nuclear  isotope  present  in  the  cell,  and  i  refers  to  the  condi¬ 
tion  that  either  a  left  or  a  right  hand  circularly  polarized  beam 
is  used  in  the  experiment.  1 s  are  the  transverse  relaxation 
times  and  p's  are  the  initial  phases.  Data  are  then  analyzed 
utilizing  a  Kalman  filter  technique  to  obtain  the  nuclear  pre¬ 
cession  frequencies. 

The  samples  used  in  the  measurements  are  all  1-ml  spherical  cells 
made  of  pyrex  glass.  The  glass  envelopes  were  baked  overnight  at 
300°C  under  high  vacuum.  After  the  bake,  rubidium  metal  was 
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129 

distilled  in;o  the  cell,  and,  a  pair  of  noble  gas  isotopes,  Xe 
83  129  13 1 

and  Kr  ,  or  Xe  and  Xe  ,  together  with  buffer  gases  such  as 
and/or  He,  were  filled  into  the  cell.  In  order  to  minimize 
the  uncertainties  of  the  gas  pressures  in  the  sealed-off  cells, 
each  cf  the  spherical  cells  was  completely  covered  with  thermal 
insulation  while  the  cell  stems  were  sealed  with  a  working  flame. 
The  rubidium  metal  in  the  sealed  cells  was  then  carefully  driven 
into  the  small  region  of  the  cell  stem  area.  The  cells  so  made 
were  then  "aged"  in  an  oven  at  80°C  for  no  less  than  two  weeks 
before  the  measurements. 


Experimental  Results  and  Discussions 


The  results  are  presented  for  the  three  different  sample  cells 

12  9 

listed  in  Table  I.  Cells  1-44  and  1-48  are  for  the  Xe  -Kr 

nq  I  qi 

oair  and  cell  1-75  is  for  the  Xe  -Xe  pair. 


129 

In  figure  8,  the  magnetic  fields,  AH,  experienced  by  Xe  and 
83 

Kr  in  cells  1-44  and  1-48  due  to  the  spin  exchange  interaction 
with  the  optically  oriented  rubidium  vapor  are  shown  as  a  function 
of  the  cell  temperature.  The  functional  form  of  temperature 
dependence  of  the  magnetic  field  qualitatively  confirms  the 
theory  discussed  earlier.  As  the  cell  temperature  increases,  the 
rubidium  number  density,  NA,  increases,  and,  at  the  same  time, 
the  rubidium  polarization,  <Sz> ,  decreases,  with  AH  being  a  func¬ 
tion  of  the  product  of  the  two,  N.<S  >.  The  magnitude  of  AH 

129  a  z 

experienced  by  Xe  is  larger  than  that  by  Kr  at  all  tempera¬ 
tures  in  the  same  cell  (i.e.  the  same  rubidium  polarization 


environment).  This  is  the  result  of  the  different  magnitudes  of 


the  enhancement  factors  (i.e. 
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Figure  8.  The  magnetic  fields  experienced  by  Xe^"2^ 
(solid  square)  and  Kr®3  (empty  square) 
nuclei  in  two  different  sample  cells  as 
a  function  of  cell  temperature 
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Comparing  the  magnitude  of  AH  experienced  by  each  isotope  in  the 
two  different  cells,  the  ones  in  cell  1-48  are  larger  thar.  those 
in  cell  1-44.  This  is  understood  to  mean  that  higher  ^  pressure 
(100  torr)  in  cell  1-48  provides  a  higher  degree  of  rubidium 
polarization  than  that  provided  in  cell  1-44  (15  torr  of  ^ ) . 

Attention  is  directed  in  figure  8  to  the  fact  that  the  ratio 
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of  the  magnetic  field  experienced  by  Xe  to  that  by  Kr  due 
to  the  spin  exchange  is  insensitive  to  the  cell  temperatures. 

The  same  ratio  is  observed  within  the  experimental  uncertainty 
even  in  the  two  different  cells.  This  is  an  expected  result  from 
Eq.  (64).  The  ratio  determined  by  this  measurement  is  2.2  ±  0.1. 

Essentially  the  same  features  are  observed  for  cell  1-75,  which 
129  131 

is  for  the  Xe  and  Xe  pair,  as  shown  in  figure  9.  The  mag¬ 
netic  fields,  AH,  experienced  by  the  two  xenon  isotopes  are 
nearly  identical  at  all  temperatures.  This  confirms  experimen¬ 
tally  the  theoretical  argument  presented  earlier  that,  under  the 
same  alkali  polarization  environment,  a  pair  of  nuclear  isotopes 
of  the  same  noble  gas  would  see  nearly  equal  magnitudes  of  mag¬ 
netic  fields  due  to  spin  exchange  with  the  optically  oriented 
alkali  vapor.  The  degree  of  the  alkali  polarization  in  cell  1-75 

may  be  close  to  that  in  cell  1-44,  in  view  of  the  fact  that  the 
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magnitude  of  the  magnetic  fields  experienced  by  Xe  is  nearly 
the  same  in  the  two  cells. 

The  same  data  shown  in  figures  8  and  9  are  presented  in  figures  10 
and  11,  respectively,  in  terms  of  NMR  gyro  bias  shift  as  a  func- 

83  i  o  q 

tion  of  temperature.  Figure  10  shows  Kr  -Xe  cell  bias  shifts 
of  -60  deg/hr  and  -100  deg/hr  for  a  1°C  variation  of  the  cell 
temperature  at  a  nominal  cell  operating  temperature  of  75°C 
(temperature  sensitivity). 
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Figure  9.  The  magnetic  fields  experienced  by  both  Xe 

and  Xel31  nuclei  as  a  function  of  the  sample 
cell  temperature.  The  difference  of  the 
fields  experienced  by  the  two  nuclear  species 
are  too  small  to  be  clearly  seen  in  the  plot. 
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Fiqure  10. 


NMR  gyro  bias  versus  the  cell  temperature  for 
the  cells  with  the  Xe^29  and  Kr83  pair  of 
isotopes.  The  same  data  shown  in  figure  3-3 
are  used. 
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NMR  gyro  bias  versus  the  cell  temperature  for 
the  cell  (1-75)  with  the  Xe-^-29  ancj  Xe^l  pair 
of  isotopes.  The  same  data  shown  in  fig¬ 
ure  3-4  are  used.  The  large  error  bars  are 
due  to  the  computational  error. 
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As  shown  in  figure  11,  the  gyro  bias  shift  in  the  cell  with  the 

Xe^^  and  Xe^^  pair  (cell  1-7  5)  is  at  least  an  order  of  magni- 
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tude  smaller  than  that  in  the  cells  with  the  Xe  and  Kr  pair. 
The  large  error  bars  in  figure  11  indicate  the  computation 
error.  (Two  nearly  equal  numbers  have  to  be  subtracted  to  get 
the  gyro  bias. )  The  existence  of  any  gyro  bias  shift  in  this 
cell  is  not  surprising  in  view  of  the  discussions  mentioned  pre¬ 
viously.  (A  few  percent  of  this  bias  shift  could  have  been 
caused  by  the  different  polarizabilities  of  the  two  circular 
polarizers  used  in  the  experiment.)  Notice  that,  at  78°C,  the 
bias  becomes  maximum,  i.e.,  when  the  cell  is  operated  at  this 
temperature,  sensitivity  of  the  gyro  bias  will  be  identically 
zero. 

3.4  SUMMARY  AND  FUTURE  PLANS 

A  theory  of  nuclear  precession  frequency  shift  due  to  the  spin 
exchange  interaction  with  an  optically  oriented  vapor  was 
reviewed.  Experimental  results  are  qualitatively  in  good  agree¬ 
ment  with  the  theory.  In  particular,  it  has  been  confirmed 
experimentally  that  a  pair  of  nuclear  isotopes  of  the  same  noble 
gas  atoms  experience  nearly  equal  frequency  shifts  under  an 
identical  alkali  polarization  environment. 
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The  gyro  bias  shift  in  the  cell  with  the  Xe  and  Xe  pair  is 

at  least  an  order  of  magnitude  smaller  than  that  of  the  cell  with 
129  83 

the  Xe  and  Kr  pair.  It  has  been  found  in  this  work  that 
there  is  a  certain  cell  temperature,  at  which  the  temperature 
sensitivity  of  the  gyro  bias  shift  is  identically  zero.  In 
order  to  possibly  utilize  the  point  of  zero  temperature  sensi¬ 
tivity  in  an  NMR  gyro,  the  existence  of  any  gyro  bias  shift 


3-23 


404577 


would  have  to  be  investigated  further  within  the  frame  of  the 
responsible  physical  phenomena  discussed  in  the  theoretical  sec¬ 
tion.  Some  of  the  possible  causes  could  be  studied  advanta¬ 
geously  by  varying  the  temperature  of  the  cell  wall  while  main¬ 
taining  the  product  of  alkali  number  density  and  the  alkali 
polarization  constant.  This  could  be  done  in  an  experimental 
arrangement  in  which  the  wall  and  the  cell  tip  (where  the  alkali 
metal  reservoir  is  confined)  temperatures  are  controlled 
separately. 
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